Resonant modes in triangular dielectric cavities 
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The resonant optical modes ol a high permittivity dielectric prism with an equilateral triangular 
cross section are discussed. Eigenmode solutions of the scalar Helmholtz equation with Dirichlet 
boundary conditions, appropriate to a conducting boundary, are applied for this purpose. The par- 
ticular plane wave components present in these modes are analyzed for their total internal reflection 
behavior and implied mode confinement when the conducting boundary is replaced by a sharp di- 
electric mismatch. Improvement in TIR conflnement by adjusting the longitudinal wavevector 
is also discussed. For two-dimensional electromagnetic solutions (kz = 0), TE polarization leads to 
longer lifetime than TM polarization, assuming that escape of evanescent boundary waves at the 
corners is the primary decay process. 

PACS numbers: 41.20.-q, 42.25.-p, 42.25.Gy, 42.60.-v, 42.60.Da 
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I. INTRODUCTION 

Micro-sized dielectric cavities are becoming increas- 
ingly important, due to their applications especially 
in microlasers and resonators with various geometries, 
including disks, ||l| triangles, Q squares, jsj |^ and 
hexagons. H, |^ Materials where the two-dimensional 
cross-section controls the properties include, for exam- 
ple, various cleaved semiconductor structures Q as well 
as zeolite ALPO4-5 crystals.]^ Even more interesting 
are lasing semiconductor pyramids that grow by natu- 
ral processes. |8[ ||, |lO| 

An understanding of the geometry dependence of the 
total internal reflection (TIR) that produces resonant 
modes can lead to improvements in designs of optical 
devices. For the reduced two-dimensional electrodynam- 
ics (fields independent of a z-coordinate along an axis 
of symmetry), the boundary element method |ll| is very 
powerful for determining the scattering fields in TM or 
TE polarizations. Alternatively, here we present a sim- 
pler analysis of the resonant modes based only on an 
analysis of the conditions needed for TIR within a dielec- 
tric cavity. The calculation is approximate but simple, 
especially for a geometry where the related Helmholtz 
equation has an exact analytic solution. 

Here we consider the modes in equilateral triangular- 
based prisms and their two-dimensional (2D) analogs; the 
triangular system is chosen here for its interesting sym- 
metries and known analytic Helmholtz solution.^, 
m, H |l|, As a first approximation, a 2D dielectric 
system with Dirichlet boundary conditions (DBC) is con- 
sidered, taking the fields equal to zero at the boundary, 
equivalent to a metallic or conducting boundary. The 
exact analytic solutions for this 2D triangular system are 



reviewed; each mode is a superposition of only six plane- 
wave components. We describe how the mode informa- 
tion is useful for determining which modes will be con- 
fined if the dielectric is surrounded by a lower index di- 
electric (or vacuum), rather than a perfect conductor, us- 
ing the conditions for TIR. A careful analysis shows that 
when a mode is strongly confined in the cavity by TIR, 
Dirichlet boundary conditions apply (approximately) to 
the relevant Helmholtz equation for both the TM and 
TE polarizations. The simple plane wave components of 
each mode can be analyzed in terms of Maxwell's equa- 
tions and the related Fresnel amplitude ratios, which are 
slightly different for TM and TE polarizations. ||l|| These 
differences are shown to imply longer estimated lifetimes 
for the TE modes, enhanced by a factor of the order of 
the squared index ratio (n/n')^, when compared to the 
lifetime of the corresponding TM mode. 

Obviously some error is made compared to employing 
a more technically complete solution of Maxwell's equa- 
tions such as boundary element methods; ||ll| the fields 
discussed here for the dielectric mismatch boundary are 
not close to the true fields unless the mismatch is very 
large. This approach, however, should approximately de- 
termine the symmetries of the confined modes and give 
reasonable estimates for the dielectric mismatch needed 
for their confinements. The improvement of confinement 
for finite-height prisms with the same cross sections, al- 
lowing for nonzero longitudinal wavevector kz , is also dis- 
cussed. 



II. SIMPLIFIED QUASI-TWO-DIMENSIONAL 
HELMHOLTZ PROBLEMS 
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Within a cavity with electric permittivity e and mag- 
netic permeability /x, we consider the solutions of a scalar 
wave equation for any component ■0 of the electric or 
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magnetic field, 

v^*-f£» = o, (1) 

where c is the speed of light in vacuum and t is time. 
The cavity shape being considered is a prism with an 
equilateral triangular cross-section of edge a in the xy- 
plane; the height is h along the z axis. Ultimately, we 
want to discuss which modes will be confined in the cavity 
of index of refraction n = ^/e/^ when it is surrounded by 
a uniform medium of lower index of refraction n' = \/e' fJ.' 
(for example, vacuum). 

We first analyze the modes, assuming the fields go to 
zero at the cavity walls, using Dirichlct boundary condi- 
tions (DBC), which corresponds to perfectly conducting 
walls, or a mirrored cavity. The range of applicability of 
these boundary conditions for a mode that is confined by 
TIR caused by an index mismatch, rather than reflect- 
ing boundaries, is discussed subsequently. A solution at 
frequency is sought, with e~*^* time dependence, cor- 
responding to three-dimensional (3D) wavevectors with 
magnitude 

K = ^, c=^; (2) 

the speed of light in the cavity medium is denoted as c* . 
Then we are solving an eigenvalue problem (Helmholtz 
equation) within the given geometry, 

VV = ~K^^. (3) 

In the general case, we can consider that we are look- 
ing for a superposition of plane waves in the form 
gi(if r-nt) ^ all of which satisfy (^, and with the cor- 
rect linear combination to satisfy the required boundary 
conditions. These different components have the same 
squared wavevector, although they can possess different 
directions of propagation and different amplitudes. The 
symmetry of the situation, however, requires only some 
reduced set of possible wavevectors for triangular-based 
prisms. 

A wavevector is composed from a z-component, k^, 
along the vertical or longitudinal axis, and the remaining 
2D components in the xy plane, = k'^ + ky, so that we 
decompose the squared magnitude as 

7^2 ^kl + k'^. (4) 

Then the frequency is partitioned accordingly, 

n^ = ujl + Lu^, (5) 

where the xy frequency lo and wavevector k are related 
by an equation analogous to (j^), as fc = w/c*. The re- 
maining part, LOz, is determined for a prism of height h, 
assuming some ip oc sin^k^z) longitudinal dependence, 
by 

kz = j, 1 = 1,2,3... (6) 



which automatically enforces DBC at the ends of the 
prism. As a result, for the prism geometry, we only solve 
the separated Helmholtz equation for two dimensions, 

V^^V = -k'i^- (7) 

A. 2D Electromagnetics at a boundary between 
media 

For 2D E&M problems, which are defined by no depen- 
dence of the fields along the prism axis, one takes kz = 0, 
ignoring any boundary conditions at the prism ends 
(which can be considered at infinity). Then Maxwell's 
equations and their associated physical boundary condi- 
tions show that one need only consider transverse mag- 
netic (TM) and transverse electric (TE) polarizations of 
the fields. (Presence of nonzero kz for a dielectric waveg- 
uide surrounded by a different dielectric medium, in gen- 
eral, does not lead to separated TM and TE modes, see 
Ref. ^.) It is well-known that in a conducting waveguide 
or resonator, at the boundaries of the 2D cross-section, 
the field ij} = Ez for TM modes satisfies DBC, and the 
field ij} = Bz for TE modes satisfies Neumann bound- 
ary conditions (NBC). Alternatively, if the medium is 
surrounded simply by a different (nonconducting) dielec- 
tric medium, there still are TM or TE polarizations, but, 
strictly speaking, these would satisfy the more involved 
boundary conditions of Maxwell's equations, rather than 
the oversimplified DBC or NBC. If the fields are strongly 
undergoing TIR (incident angle well beyond the critical 
angle), however, both polarizations are shown here to 
satisfy DBC, to a certain limited extent. This can be 
seen by examining the actual vector fields present under 
the refiection and refraction at a boundary between two 
media, for the two possible polarizations. 

Consider a planar boundary between two media, where 
the boundary defines the xz-plane. The region y < is 
occupied by a medium of index n = ^/fjl, while the re- 
gion y > is occupied by a medium of index n' = ^t' pJ , 
with n > n'. Primes refer to quantities on the refracted 
wave side (ultimately, outside the system we are study- 
ing). The wavevector magnitudes are k = ^V^M ^-^d 
k' = ^ y/ e'/z' in the two media. All the waves in this 
discussion are assumed to have e~*'^* time dependence. 
A 2D plane wave '-^ ^^ik^^+kyy) t^^Ii fields Ei, Bi, propa- 
gating in medium n with ki = {kx,ky) = k[sui9i, cos 9 i) 
and incident on the boundary, can be polarized in two 
ways, which correspond exactly to the TM and TE po- 
larizations. 

Diagrams of the two possible situations are given, for 
example, in Chap. 7 of Ref. |l8[ For Ei perpendicular 
to the plane of incidence {xy-plane), Ei = EiZ, and the 
magnetic field will have no z-component anywhere, which 
is the same as the TM polarization. On the other hand, 
for Ei lying within the plane of incidence, Bi = BiZ, and 
the electric field has no z-component anywhere, corre- 
sponding to the TE polarization. Thus, we can use the 
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Fresnel reflection and refraction amplitudes for these two 
cases, that result from the actual boundary conditions for 
Maxwell's equations, to investigate the appropriateness 
of applying DBC for fields undergoing TIR. Furthermore, 
this analysis will be used for determining the properties 
of the evanescent refracted wave under TIR; the power 
lost in this boundary wave can be expected to play a role 
in the lifetime of a resonant mode. 

TM polarization: For incident Ei = E'OfgjCs^a^+fc^y) 
perpendicular to the plane of incidence, there is also a 
reflected wave Er = E'^ze^'^^''^~^^y^ with the same polar- 
ization, but a different magnitude and phase, depending 
on the angle of incidence 6*^. One can express the reflected 
wave amplitude E^ by the Fresnel formula: 



= E^ e- 



(8a) 



threshold for TIR [6 = 6*^), (|Tl|) gives a = 0, whereby 
( |l^ ) indicates that the E^ field now will peak at y = 
and satisfies a Neumann BC. 

Seeing that the BC on E-^ ranges between the two ex- 
tremes of DBC and NBC, clearly neither boundary con- 
dition is fully applicable to the TIR regime. However, 
it suggests that one should use NBC for finding the TIR 
threshold conditions, and DBC to determine the field dis- 
tributions at large enough index mismatch. One can get 
some indication of the crossover to DBC-like behavior by 
locating the incident angle 6i = 9 where the phase shift 
passes ^, i.e., when tan ^ = 1. From (11) one finds 



cos 9 



cosf 



(14) 




Here 9' is the angle of the refracted wave, obtained from 
Snell's Law, n sin 9i — n'sinfl', which implies TIR when 
9i surpasses the critical angle 9c, defined by 

n' 

sin6'c = — . (9) 

Under TIR, E° has the same magnitude as E^, but 
is phase shifted by angle a, because the cosine of the 
refracted wave becomes pure imaginary: 

cos6'' = i7, 7 = \/(sin6li/ sin 61^)2 - 1. (10) 

Then it is useful to express the phase difference between 
the incident and reflected waves as 

tan^ = ii./22£!^ (11) 

2 ^' Y C0S2 9^ ^ ' 

The linear combination of incident and reflected waves 
in medium n is E ^ Ei + Er ~ EzZ, having the spatial 
variation approaching the boundary (region y < 0), 

E,=E, + Er^ 2i;^-''te**'--" cos (fc^y + |) (12) 

Corresponding to this is the associated evanescent wave 
in medium n' (region y > 0), 

E', = 2E'^e-''^e'^-''e~^'''y, (13) 

where = k'^ due to Snell's Law. Inspection of Eq. 
( p^ ) shows that the E^ held on the incident side acquires 
a node at the boundary y — Q and satisfies Dirichlet 
BC only when the phase shift attains the value a = tt. 
According to (pl|), this occurs only in the limit 9i 
90°, i.e., extreme grazing incidence. Alternatively, at the 



In the usual practical situation, with /i = /i' w 1, we have 
cos ^ = cos 9c- For usual optical media with very weak 
magnetic pro perties, one sees that a large index ratio 
n/n' m does not increase the range of applicability 

of Dirichlet BC for TM polarization. To give a numerical 
example, for a weak index mismatch with sin^^ = 1/2, 
giving 9c = 30°, the crossover angle is 6* = 52.2°; in 
order to reach phase angle a = 0.9tt, quite close to DBC, 
requires an incident angle di — 82.2°. Even for a larger 
mismatch sin6'c — 1/4, with 6c — 14.5°, one gets only a 
slight improvement to 9 — 46.8°, and to get to a = 0.97r 
still requires 9i =81.3°. 

The conclusion is that it makes some reasonable sense 
to apply NBC to find the limiting conditions for TIR, 
but, in general, provided the incident angle is sufficiently 
larger than the crossover angle 9, an approximate descrip- 
tion of the fields on the incident side should be possible 
using DBC. The application of DBC to this problem im- 
proves with higher index mismatch, but not as strongly 
as one would hope, because the phase angle does not de- 
pend on the dielectric permittivities for this polarization. 

TE polarization: The discussion is similar, but now 
the magnetic field B is polarized in the z direction ev- 
erywhere, and controls the other fields, according to re- 
lations B = y/ejl k X E tor each plane wave, and am- 
plitude relations B = ^/em Taking incident wave 
B, = 50£gj(fe,a;+fc„y)^ ^-^^-^ elcctric field amplitude E° = 
B'j'/y/eJI, there is a refiected wave Br = i?°ze'(''"^~''!'^', 
with electric field amplitude E^ = B^/ ^/ejl. Now the am- 
plitude ratio E^/Ef — e~*" is described by the different 
Fresnel formula: 




Clearly, the same amplitude ratio also applies to B^ j B^ . 
The phase difference between the incident and reflected 
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waves can be expressed as 



a e cos^ Oc 



1. 



(16) 



The linear combination of incident and reflected magnetic 
waves in medium n is B — + Bj. — B^z^ and mirrors 
the behavior of E for the TM problem, having the spatial 
variation approaching the boundary, 



Bz — Bi + Br 



2B^e 



^cos(fcj,y+|) (17) 



The associated evanescent wave in medium n' is 



B' = 2B°e' 



t 2 ^ikxX ^• 



(18) 



Clearly, the behavior of Bz near the boundary for TE po- 
larization is the same as that for near the boundary 
for TM polarization. It means that provided the incident 
angle is far enough beyond the critical angle, one should 
also apply DBC for TE fields; NBC would only be rea- 
sonable just beyond the TIR threshold. However, due 
to the presence of the permittivity ratio in (|T^), a large 
index mismatch does enhance the applicability of DBC 
for TE polarization. The crossover incident angle 9 at 
1 is now 



which tan ^ 



cos 9 — 



cost 



cos 9r 



1 



(19) 



where the latter expression applies when fi ^ fi' . For 
some numerical examples, now the modest index mis- 
match with sin 0c = 1/2 and 9c = 30° leads to a very 
nearby crossover angle 9 = 32.8°, meaning that the 
strong TIR regime and region of adequate applicability 
of DBC is very wide. For larger mismatch sin 0c = 1/4 
with 9c — 14.5°, the DBC regime is even closer to 9c, be- 
ginning around 9 = 14.9°. Thus, application of Dirichlet 
BC to finding the 2D resonant modes for a cavity with 
TE field polarization should be very acceptable, even 
more so than for TM polarization, except when the mode 
has plane wave components extremely close to their TIR 
threshold. 

Within the limitations indicated in the above analysis, 
we continue by discussing the modes in a 2D equilateral 
triangle, under the assumption of DBC for either the TM 
or TE polarizations. 



Exact modes for an equilateral triangle with 
DBC 



The triangular cross section offers an opportunity for 
exact solutions to the 2D Helmholtz equation. This 
problem has been solved analytically p^, f3l f4l 
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for 

both DBC and NBC in several contexts, including quan- 
tum billiards problems, [p7| ITsI, |20|, pT| quantum dots, pi] 




FIG. 1: Description of the coordinates for a 2D triangular 
cavity of edge a. The geometrical center at C, the origin of 
the a;i/-coordinates, is a distance above the lower (6o) 
edge. The lower left corner at L is the origin of the skew XY- 
coordinates. The dashed line demonstrates the reflections of 
a ray originating at angle (3 = 70° to the lower edge, requiring 
two complete circuits to return to the same angle. 



and lasing modes in resonators and mirrored dielectric 
cavities. ^ In particular, it is interesting to note that 
only the triangles with angle sets 7r(i,i i), 7r(i,|;,;|), 
and 7r(^, I, g) are classically integrable[pO| and have sim- 
ple closed form wavefunction solutions derived in various 
ways,[p^ 14, 21 following the first solution for trian- 
gular elastic membranes by Lamc.|ia Here we use the 



equilateral triangle of edge length a for its interesting 
symmetries and resulting simplifications. 

Coordinates are used where the origin is placed at the 
geometrical center of the triangle, and the lower edge is 
parallel to the x-axis, as shown in Fig. |^. The notation 
6o, ^1, and hi is used to denote the lower, upper right, 
and upper left boundaries, respectively. This discussion 
concerns DBC. 

Following Chang et aL,|^] some comments can be made 
on the sequence of reflections a plane wave trapped in the 
cavity makes, but with a slightly different physical inter- 
pretation. A plane wave leaving the lower face (6o) at 
angle /3 to the x-axis sequentially undergoes reflections 
at the other two faces, generating plane waves at angles 
240° - (3, /3- 120°, relative to the x-axis. Finally, when 
reflected again off the lower face, it comes out at 120° — /?, 
which does not match the original wave unless (3 = 60°. 
However, when allowed to propagate again around the 
triangle, the sequence of angles is 120° -I- /?, — /3, which 
then comes out at +(3 after the reflection off the lower 
face. So the wave closes on itself after two full revolu- 
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Lowest Modes 

2D Triangular Cavity 



2 4 6 8 

m 

FIG. 2: Scaled frequencies of some of the lowest possible 
modes in a 2D triangular system, as a function of quantum 
index m, for indicated n > m. Indices n and m must be 
integers of equal parity, c* is the light speed in the medium. 



where the rotated vectors are 



Rki 



R% 



Rk2 



R^k2 




X 11 

2 2\ 




(23) 



(24) 



By design, these rotated waves satisfy the boundary con- 
ditions, ■01 (^i) = and i/'2(fe2) = 0. 

In order to determine the allowed (fci,fc2), one needs 
to impose DBC on all three boundaries, for a linear com- 
bination with unknown coefficients ^0,^1,-^2, 



1/' = Aoipo + Aiipi + A21P2 



(25) 



tions, no matter what the initial angle. As an example, 
the sequence of reflections starting with P = 70° from 
boundary bo are shown in Fig. |l]. In general, any origi- 
nal wave simply generates a set of six symmetry related 
waves rotated by ±120° and inversions through the y- 
axis. 

These considerations show that the general solution is 
a superposition of six plane waves, which can be obtained 
by 120-degree rotations of one partially standing wave. 
Consider initially a wavevector with xy Cartesian com- 
ponents, k = (fci, fc2), and defined basic vectors ki = kix, 
k2 — k2y, producing a wavefunction written as 



V'o 



k2a 
2V3 



(20) 



This is a combination of plane waves at two angles /3 
and —(3 to the x-axis, where tan/? = ^2/^1, and the 
combination goes to zero on boundary 60, where y = 
~2^' ^^^^ ^^^^ ^ standing wave. By assumption, 
we need ^2 7^ to have a nonzero wavefunction. As 
discussed above, when refiected off the triangle faces bi 
and &2, this simply produces rotations of ±120°. If R 
is an operator that rotates vectors through +120°, the 
additional waves being generated are 



iPi = e^^-^'^i^^^sin 



sm 



iRk2 



{R^k2) 



k2a 



2V3. 



k2a 



2\/3. 



(21) 



(22) 



Imposing DBC on all boundaries determines the allowed 
wavevector components as 



ki = —m, 
6a 



da 



m = 0, 1,2.. 



1,2,3. 



(26) 



(27) 



iTrm _ gjTrn appears; 



Furthermore, the parity constraint e 
that is, n and m are either both odd or both even. 
The resulting xy frequencies are given by 

c 27r 



kl 



e/i 3a 



\/ m? + 3n^ 



(28) 



The mode wavefunctions are described completely using 
the amplitude relationships that result: 



A2 = Ane 



(29) 



Therefore, solutions are specified by a choice of integers 
m and n and the phase of the complex constant ^q. In 
Fig. H the frequencies of the lowest modes are presented, 
scaled with the speed of light in the medium, c* , and ar- 
ranged as families for each value of the quantum number 
n. This solution is represented in the physically moti- 
vated form described by Chang et a/.,||] and is entirely 
equivalent to the first solution given by Lame'|jl^ and 
revisited by various authors, [p^ 

The solutions obtained have obvious symmetry prop- 
erties. We can get all the possible eigenfrequencies by 
applying the restrictions, < m < n. Choices with 
TO > n or m < also give allowed frequencies, how- 
ever, these only correspond to other k rotated by ±120° 
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(b) mode 4 



2nd excited state 
(nondegenerate) 




m=0, n=4 
coa/c*=14.5104 
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FIG. 3: Wavefunctions of the two lowest modes with m = 
in a 2D triangular system: the ground state (a), with n = 2, 
and the third excited state (b), with n = 4, at double the fre- 
quency of the ground state. Solid/partial shadings represent 
positive/negative wavefunction values, whose magnitudes cor- 
respond to the radii of the symbols. A grid with N = 40 was 
only used to present the diagrams. Grid sites without symbols 
have \tp\ < 0.021V' max - 



from some original k defined with < m < n. For in- 
stance, having found k = (^1,^2) = |^(m, VS^i), values 
of k' = (fc'i, fc^) = |^(m', \/3n') corresponding to ±120° 
rotations result from 



k[ = fcicosl20°Tfc2sinl20° 
k'2 = ±ki sin 120° + k2 cos 120° 



(30) 



FIG. 4: Wavefunctions of the doubly degenerate first excited 
state in a 2D triangular system. Two orthogonal sub-states 
corresponding to different choices of the phase of ijj are pre- 
sented in (a) and (b). 



which implies transformation to the new mode indexes, 



m 



1 3 

2 2 



1 

±-m 



(31) 



As specific examples, the choice (m, n) = (1, 3) gives one 
mode, which after rotations by ±120° could also be de- 
scribed by the integers (m, n) = (4,2) or by (m,n) = 
(5,1). Note, however, this represents only one mode, 
with different choices of the reference edge or triangle 
base. Wc should also note, any mode with m 7^ is 
doubly degenerate; the change m —m gives an inde- 
pendent mode with the same frequency, which rotates in 
the opposite sense around the triangle. A real represen- 
tation of the degenerate pairs comes from taking the real 
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(a) mode 5 

3'''^ excited state / \ m=2, n=4 
1st wavef unction / • •/,•.•. \ Q)a/c*=15.1029 
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The wavefunctions of the lowest modes with m = 
are displayed in Fig. ||. All the modes with m — have 
this typical triangular shape, periodically patterned over 
the entire system. Wavefunctions of the first two ex- 
cited states are presented in Figs. ^ and |^. Each of these 
states, having nonzero m, is doubly degenerate, therefore, 
there is some arbitrariness in the wavefunction pictures. 
Two degenerate wavefunctions were obtained from the 
real and imaginary parts of ■0 in Eq. (^); an equivalent 
degenerate pair can be obtained by changing the complex 
constant Aq. The degenerate pairs can be shown to be 
orthogonal in the usual sense. Wavefunctions for the se- 
quence of many of the other lo west modes are displayed at 
www.phys.ksu.edu/~wysin/. The modes described here 
are the complete set of allowed modes for the triangular 
geometry. 



(b) mode 6 

3'''^ excited state 
2nd wavefunction 
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C. Equilateral triangle with NBC 

Although we do not apply these in this work, it is inter- 
esting to note the minor differences in the solutions when 
the triangle is solved with Neumann BC. As far as the 
changes in the wavefunctions, NBC simply requires the 
sine functions in the expressions for -00, ^i, and ■02, Eqs. 
@, (pi), to be replaced by cosine functions. The 

expressions ( |26| ) and for the allowed wavevectors are 
still valid, except that the restriction m < n should be 
changed to m < n. So the spectrum is expanded slightly, 
there is an entire set of nondegenerate lower modes with 
m = n (alternatively with rotated indexes (m, 0)); the 
NBC ground state is lower than the DBC ground state. 



D. 3D prisms 



FIG. 5: Wavefunctions of the doubly degenerate third excited 
state in a 2D triangular system. Two orthogonal sub-states 
corresponding to different choices of the phase of are pre- 
sented in (a) and (b). 



and imaginary parts of ip to form two independent wave- 
functions. Different values of the complex constant 
will produce other choices of the two independent wave- 
functions. For m — 0, the wavefunction can be made 
pure real; the modes are non-degenerate and n must be 
even, since m and n must be of equal parity. 

The ground state (Fig. ||a) has (to, n) — (0, 2) and cor- 
responding xy frequency, cu = -^|- . This corresponds 

to a wavelength of A = ^ = , which is equal to the 
height of the triangle. Naively, one might have expected 
that a half-wavelength of the lowest mode could have fit 
across the triangle, as would occur in simple rectangular 
or one-dimensional geometry. Apparently, that choice 
would not have the ability to fully satisfy DBC on all 
three surfaces. 



For a 3D triangular-based prism of height h, and base 
edge a, the solutions to the wave equation can be written 
as products of an a;?/-dependent part ipxy and a part ipz 
depending only on the vertical coordinate z. The cor- 
responding k and kz have already been discussed above. 
Solutions can be found combining k and ipxy from the 
2D solutions with the k^ and ij^z discussed above. This 
approach would assume DBC over the all surfaces of the 
cavity. While the physical situation would not exactly 
satisfy the requirements for classification of modes in TM 
and TE polarizations, the error in this assumption should 
be least at small k^. Thus, it could give some rough idea 
about the modification of the confined modes due to the 
finite vertical height of the cavity. 

III. TIR CONFINEMENT AND THE 
RESONANT MODE SPECTRUM 

Knowing the eigenmodes of the cavity, various ques- 
tions can be addressed and answered based on the struc- 
ture of the wavefunctions. The primary question is: 
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which modes are confined by TIR (or resonant) when the 
previously assumed Dirichlet boundary conditions are re- 
placed by an index mismatch at the boundary? As seen 
below, some states can never be TIR confined. States 
which have the possibility for confinement will be refered 
to as TIR states. Related to this is to estimate a decay 
time for any of the TIR modes. 

As suggested in the Introduction, the mode confine- 
ment can be determined approximately by identifying 
those modes which completely satisfy the elementary re- 
quirements for total internal reflection. We assume that 
outside the system boundary a different refractive ma- 
terial with relative permeability /x' and permittivity e' 
(which could be vacuum for greatest simplicity) and in- 



dex n' = 



rather than a conducting boundary. The 



wavefunction analysis above describes a mode as a su- 
perposition of plane waves with the specified (full 3D) 
wavevectors K. Each plane wave component can be an- 
alyzed to see whether it undergoes TIR at all system 
boundaries on which it is incident. If all the plane wave 
components of a mode satisfy the TIR requirements, then 
the mode is confined; it should correspond to a resonance 
of the optical cavity. If any of the plane wave components 
do not satisfy TIR, the fields of the mode will quickly leak 
out of the cavity and there should be no resonance at that 
mode's frequency (also, without TIR, the assumption of 
DEC is completely invalidated) . These results clearly de- 
pend on the index ratio from inside to outside the cavity, 
denoted as 



(32) 



A plane wave within the cavity with (3D) wavevector 
Ki has an incident angle on one of the boundaries ex- 
pressed as 



sin 9i 



K,, 



(33) 



where || indicates the component parallel to the bound- 
ary. TIR will take place for this wave provided 



sm Oi > sm ( 



1 

N' 



(34) 



For 2D geometry, Eqs. (^) is applied directly to 2D 
wavevectors, Ki = ki: one just needs to identify the 
parallel component of fcj relative to a boundary, for each 
for each plane wave present in ■0, Eq. (^). The main 
difhculty is to satisfy TIR on all boundaries on which 
that component is incident. 

For 3D prismatic geometry, the wavevector Ki includes 
a ^-component. The net parallel component needed in 
( ^ ) can include planar (xy) and vertical (z) contribu- 
tions. The addition of a ^-component will have the 
tendency to raise the energy of any mode, and the net 
wavevector magnitude, thereby improving the possibility 
for mode confinement. 



A. Confinement in 2D triangular systems 

For both the TM and TE 2D modes of an equ ilate ral 
triangle, we can use the exact solutions in Sec. II B to 
investigate whether these can be confined by TIR. Due 
to the threefold rotational symmetry through angles of 
0°, 120° and -120°, the TIR condition need only be ap- 
plied on one of the boundaries, say, the lower boundary 
bo, parallel to the a;-axis. The first wave ipo is com- 
posed of two traveling waves, one incident on bo and 
one reflected from bo', both have wavevectors whose x- 



component is fci 



3a 



m, where to = 0, 1, 2, .... The other 



rotated waves ipi ^-nd ?/;2 are composed from traveling 
waves whose wavevectors have x-components of magni- 
tudes i(\/3fc2±A:i), where ^2 = §f V^'t- with n — 1,2, 3.... 
For the allowed modes, to < n, or fci < /c2, which shows 
that 



ki <^{V3k2- h) <^{V3k2 + ki). (35) 

The x-component of the ipo wave is always the smallest; 
this wave has the smallest angle of incidence on bo, so if 
it undergoes TIR then so do i/ji and ■02 , and the mode is 
confined. As ki is the component parallel to boundary 
bo (i.e., fci||), the condition for TIR confinement of the 
mode is 



(36) 



or equivalently in terms of quantum numbers to and n, 



> 



N2 



1 



(37) 



This relation contains some interesting features. First, 
since all modes have m < n, the LHS is always less than 
1, and confinement of modes can only occur for ade- 
quately large refractive index ratio, N > 2. For N < 2 all 
modes will leak out of the cavity; they cannot be stably 
maintained by TIR. Secondly, for any particular value 
N > 2, relation ( p7| ) determines a critical m/n ratio; 
modes whose m/n ratio is below the critical value will 
not be stable. As m/n relates to the geometrical struc- 
ture of the mode wavefunction, there is a strong relation 
between the confined mode structures and the refractive 
index. An alternative way to look at this, is that each 
particular mode will be confined only if N is greater than 
a specific critical value determined by the m/n ratio for 
that mode: 



N > Nc 



I. 



(38) 



The relation shows that generally speaking, modes with 
smaller n/m are most easily confined; stated otherwise, 
the modes where to is closest to n are the ones most 
readily confined. On the other hand, all modes with to = 
leak out, because these have wave components with a 
vanishing angle of incidence on the boundaries. 
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FIG. 



TIR mode confinement limits (index ratio N 



\/ efi/e'fi') for 2D triangular cavities. Modes are confined 
where the m/n ratios (as indicated) lie above the solid curve, 
corresponding to TIR on all boundaries [relation (^^]. Inter- 
sections on the N-axis give the critical index ratios for each 
mode. 



Fig. ^is a kind of phase diagram for mode confinement 
by TIR. Some of the modes' m/n ratios are plotted vs. 
refractive index ratio N, together with relation (^7|). A 
particular mode is confined only if its m/n value falls 
above the critical curve. In this way one can easily see the 
critical refractive indexes for each mode. As mentioned 
above, modes where m is close to n require the smallest 
refractive index values for confinement. 

One more application of these results is shown in 
Fig. ^ where the frequency of the lowest confined mode 
is plotted versus y^e/I, with vacuum outside the cavity 
(e' fi' — 1). The quantum numbers (m,n) for the lowest 
confined mode are indicated in each curve segment. The 
two different curves correspond to plotting cua/c (dashed) 
and frequency scaled with refractive index, {uja/c)y^ejl 
(solid). Of course, as ^/ejl approaches the value 2, the 
lowest frequency becomes large and the minimizing mode 
has m very close to n, with both large. At the opposite 
limit of large values of ^/ejl, the lowest frequency becomes 
small, and the mode (1, 3) is always the lowest frequency 
mode that is confined for any -y/e/I > 5.29 . The graph 
only shows the minimum confined frequency at the inter- 
mediate values of refractive index Ve/I. 



B. Confinement in triangular based prisms 

We consider a vertical prism of height h with a triangu- 
lar base of edge a. With perfectly reflecting end mirrors 
at z = and z — h the allowed longitudinal wavevec- 
tors would be kz = /tt/Zi, where I is an integer. For a 
plane wave component with 2D wavevector fc^, the net 



FIG. 7: Frequency of the lowest confined mode for a 2D 
triangular system surrounded by vacuum, as a function of 
the refractive index. Pairs (m, n) indicate some of the modes' 
quantum numbers. No modes are confined for ^/e/I < 2. 



3D wavevector magnitude within the cavity is 



- JB + kl 



(39) 

This increases the incident angles on the walls of the 
cavity: the presence of nonzero k^ should enhance the 
possibility for confinement, compared to the purely 2D 
system. 

At the lower and upper ends of the prism, the parallel 
wavevector component of any of the plane waves present 
is the full 2D ki, which can be expressed as 

27r 
3a 

Finding the incident angle by (^ and Snell's Law (jsj), 
this implies a critical index ratio needed for TIR confine- 
ment by the cavity ends, 



— V 'Tl 



3n2 



(40) 



|\] (ends) 



1 




27r / m^ -f 3n2 



(41) 



On the vertical walls of the prism, the parallel wavevec- 
tor component is a combination of k^ and the 2D /cj 1 1 . 
If TIR occurs on on one wall then by symmetry it will 
take place on all the walls. Considering the 6o wall, the 
-00 wave has the largest incident angle as in the 2D prob- 
lem, and both ki and k^ determine its parallel wavevector 
component. 



Then this determines the critical index ratio needed for 
TIR by the cavity walls. 



'27r\' 



|\| (walls) _ 




3n2 



m^ + {4^f{kza) 



(43) 
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FIG. 8: Triangular prism critical index ratios for TIR from 
the walls (open symbols) and ends (solid symbols) of the cav- 
ity, for the modes with 2D quantum indexes m = 1, n = 3, 
and longitudinal wavevector k^,. The symbols show the al- 
lowed k-^a values Ina/h when h — a. 



For TIR mode confinement, the actual index ratio N must 

be greater than both N^''"'^'''' and N^™''*'''. 

Some results are shown in Fig. ^ for the to = 1 , n = 3 
modes as a function of k^. The symbols indicate the 
allowed kz values for the prism height equal to the base 
edge, h = a. The limiting 2D critical ratio is seen at 
the point fc^a = 0. The critical index ratio is lowered 
from the 2D value Nc ~ 5.29, down to values as low as 
Nc ~ 1.40 for kztt « 11. A more significant reduction in 
Nc occurs for the ground state, m — 0, n — 2, Fig. ||. 
It is not possible to confine this mode in 2D; whereas, it 
has Nc ~ 1.41 for k^a « 7.25. Similar significant critical 
index ratio reductions occur for all the modes tested. 

One sees that Nc*'™''*'^ and Nc^'^"'^'''* always cross at some 
intermediate k^, where they produce the smallest index 



ratio N 



(min) 



needed for TIR. Setting them equal, we find 



Nj-'") = ,/l+(l + £j 



which occurs at 



kt = 



^ ) VSn = k2. 



(44) 



(45) 



It is a rather intriguing result; choosing h/a such that a 
mode will occur at this kz will be the optimum choice for 
having the mode confined most easily by TIR. This could 
be useful for control over selection of desired modes in a 
cavity. 

Cautionary comments are in order. The above dis- 
cussion applies exactly only to a scalar wave. For elec- 
tromagnetic waves, it applies to either the TM or TE 



FIG. 9: Triangular prism critical index ratios as explained in 
Fig. ^ for the modes with 2D quantum indexes m = 0, n = 2 
(2D ground state). 



polarizations in an approximate sense for quasi-2D EM 
modes, requiring small longitudinal wavevector fc^a < 1. 
The greatest enhancements in TIR were found to occur at 
values kzQ > 1, primarily because a reasonably large kz 
is needed in order to satisfy the TIR requirements at the 
ends. The presence of significant kz values, however, will 
lead to a mixing of the TM and TE polarizations, making 
this calculation invalid. Then, for practical purposes, the 
calculation is mostly interesting for how it indicates the 
improvement in TIR confinement expected mainly on the 
2D walls of the cavity, at small kz. 



IV. TIR MODE LIFETIMES 

For the 2D solutions {kz = 0), it is interesting to esti- 
mate the lifetimes of the TIR-confined modes, contrast- 
ing the results for TM and TE polarizations. When all 
of the plane wave components in ^ satisfy the TIR con- 
ditions, there is still the possibility for the cavity fields 
to decay in time. Clearly, we have only an approximate 
solution, since DBC is not exactly the correct boundary 
condition. The effect this causes is difficult to estimate. 
Another source of decay are diffractive effects: the finite 
length of the triangle edge and the presence of sharp cor- 
ners is likely to have special influence on the TIR that is 
difficult to predict. One feature, however, which can be 
considered as due to diffraction, is the leakage of bound- 
ary waves at the corners of the triangle. [p3| Under condi- 
tions of TIR, an evanescent wave exists within the exte- 
rior medium, decaying exponentially into that medium, 
and moving parallel to the cavity surface. When it en- 
counters the corner of that edge, a sharp discontinuity in 
the surface, it can be expected to constitute power radi- 
ated from the cavity. Here we consider the mode lifetime 
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estimates based solely on the losses due to these bound- 
ary waves. 

Based on the ratio of the total energy U stored in the 
cavity fields, compared to the total power P emitted by 
the boundary waves from all the corners, an upper limit 
of the mode lifetime can be estimated as 



U 



(46) 



The calculations of U and P have slight differences for 
TM versus TE polarization. Therefore, there is no reason 
to expect these lifetimes to be the same. 

Cavity energy: For both po lariz ations, we use the 
wavefunction ip reviewed in Sec. [IB, which can be ex- 
pressed more succinctly as 



iP ^ Ao{ e^'^'^'sm k2{y 



Jki{-^x+^y+a) , 



(47) 



, la. 



The total energy of the fields within the cavity of height 
h can be written as 



U = h dxdy 



e\E\ 



h dxdy 



\Bl 
87r/i 



(48) 



the first form is convenient for TM modes (|i?P = ji/'P), 
the second is convenient for TE modes {\B\'^ = IV'P)- So 
both calculations require the normalization integral oiip. 
This integral can be simplified by a transformation to 
a skew coordinate system whose axes are aligned to two 
edges of the triangle as shown in Fig. |l]. Placing the origin 
of the new coordinates {X, Y) at the lower left corner of 
the triangle, with X increasing from to a along edge bo, 
and Y increasing from to a along edge 62 , the definition 
of the new coordinates results from 



X-Kycos60° = X + -, 
r sin 60° = 



2V3' 



(49) 



The numbers (|, j^) are simply the displacement of the 
origin (vector from triangle corner L to center C). The 
wavefunction is now expressed as 



ip = Aole 



iki{x+^Y- 



V3 



koY 



V3 



k2{-X-Y + a) 



+ e 



(50) 



In this form, it is more obvious that each term goes to 
zero on one of the boundaries, X = (62), Y = (bo), 



or X -I- F = a (&i). Using the periodicity of ip, the inte- 
gration over the triangular area is effected by 



dx dy 



dX / dY \J\ 
"'0 



(51) 



where the Jacobian is | J| — ^ and the factor of ^ cancels 
integrating over two triangles. The absolute square of V' 
involves three direct terms (squared sines involving only 
fc2) and six cross terms from Eq. (|5^). It is possible to 
show that the cross terms integrated over the triangular 
area all are zero, due to the special choices of allowed ki 
and /c2 given by ( p6[ ) and (|2^). The remaining nonzero 
parts result in 



dx dy \ipf 



3V3 



a'\Ao\ 



(52) 



Boundary wave power: The symmetry of the wave- 
function causes the boundary wave power out of each 
edge to be the same, therefore, we calculate that occur- 
ring in edge Bq (at ?; = 0) and multiply by three for the 
total power. This calculation follows that presented by 
WiersigQ for resonant fields in a regular polygon. 

Looking at the wavefunction (^^, one can see that 
there are three distinct plane waves incident on bo. First 
is the wave with the smallest angle of incidence, resulting 
from the first term in (E^, 



2i 



(53) 



Using the allowed values for ki and /c2 , the angle of inci- 
dence is seen to be 



m 



sm f n = 



(54) 



Next, there is a wave with the largest magnitude incident 
angle, due to the second term in (E^), 



Ao 



whose incident angle is 



sin6i7 



1 — m — 3n 

2 V + 3n^ 



, (55) 



(56) 



A negative value of 0i+ means the wave is propagating 
contrary to the x — axis. Finally, the last term in (jij 
leads to a wave with an intermediate incident angle, 



r2 



2i 



whose incident angle is 



sme/o = -- 



1 — m + 3n 



2 Vm-^ + 3n^ 



(57) 



(58) 



The plus/minus superscripts on these waves refer to the 
positive/negative exponents in the sine functions of Eq. 
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Now, for each of these incident waves, there is a corre- 
sponding evanescent wave propagating along the edge of 
the cavity; these are assumed to produce emitted power 
when encountering the triangle corners. The Poynting 
vector S' associated with a single plane evanescent wave 
along the bo boundary is 



S' = — 5R(-B' X H'*) = —J — \E'\''sme' x (59) 



On the other hand, the linear superposition of the three 
waves ipQ , tp^ , ipt leads to an interference pattern both 
within the cavity, and in the evanescent waves and exte- 
rior power flow. Careful consideration of a linear com- 
bination of two waves shows that, although interference 
leads to a spatially varying S" with components both par- 
allel and perpendicular to the boundary, an integral num- 
ber of wavelengths of that pattern fits along the edge. 
Thus, it is clear that the interference effects can be ig- 
nored in the calculation of the emitted boundary power. 
For the total emitted power due to boundary waves, it 
is sufficient to sum the individual powers for the three 
independent incident waves. 

For TM polarization, with -0 = i?^, the exterior 
electric field of a single eva nesce nt wave has only a z- 
component like that in Eq. ( II A ). Applying Snell's law, 
and integrating the Poynting vector from y ~ to y — oo, 
the power flow along x in one boundary wave, on one 
edge, is 



P. = 



ch \e: 



0|2 



n sin 9i 



47r^' uj/c ^(n sin 6*02 _ (n') 



2 

COS — 
2 2 



(60) 



where a is the phase shift given by (|ll|) in Sec. 



One can see that the boundary wave power has a depen- 
dence on (sin^i — sin^c)"^^^- Although each of the waves 
ipQ,ip^,ip2 ^ill produce a boundary wave emission, the 
wave has the smallest angle of incidence, and pro- 
duces by far the largest boundary wave power. This can 
be seen by examining the expressions for incident angles 
9q , of and and using the facts that m < n and m, n 
have the same parity. Therefore a good lifetime estimate 
can be made using only the power due to ipQ . From the 
expression (|5^ ) for , the squared magnitude of its elec- 
tric field is = il-^oP- From expressions ( |5^ ) and 
(^), the total TM cavity energy is 



Utm 



eh 3V3 2| . |2 



(61) 



The estimate of the lifetime due to boundary wave emis- 
sion only, from all three edges combined, is ttm ~ 
Utm/-^Px- It is convenient to express the reult in di- 
mensionless form, scaling with the mode frequency. 



WTtm 



VS /way 



1 - (sin6»c/sin6'o")2 



cos^ 



sm 



(62) 



where 9q depends on the mode quantum numbers ac- 
cording to Eq. (^J). In the usual case where ^ « the 
second line of the formula simplifies to just cos^ 9c. 

Obviously, when 9q approaches 9c, which would oc- 
cur at weak enough index mismatch, the estimated life- 
time Ttm 0, which is the limit of a non-bound state. 
At the opposite extreme of large index mismatch where 
n/n' ^ 1, this estimate varies as 1/cos^9q, and since 
Lua/c* is a number of order unity, the order of magnitude 
is determined by 



C 



(63) 



The result is interesting because it shows a lifetime that 
increases with the triangle size, as well as being propor- 
tional to the refractive index in the cavity. 

For TE polarization, with ip = B^, the exterior 
magnetic field of a single e vane scent wave has only a z- 
component like that in Eq. ( II A ) . The calculation follows 
the same reasoning as used for the TM modes, but the 
specific details lead to a slightly different result. Defini- 
tion of the Poynting vector as in ( |59| ) is the same. The 
substitution \E'\ ^ \B'J/y/7JF together with Eq. ( |IA| ) 
for B'^, followed by = within the cavity, 

leads to the boundary wave power with extra factors. 



a 



infi'e'fi uj/c y/{nsm9,)^ - (n')' 



: COS 



(64) 



where now the phase shift a given by (|l6|) in Sec. 11 A 
depends on the ratio e/e' ^ 1 rather than ^/fj,' w 1. 
Furthermore, using (52) and (48), the energy stored in 



the cavity fields is now 



Ute = 



^ ^^„2| . ,2 



(65) 



With cavity magnetic field strength |-Bf P = i|^oP, and 
again estimating the lifetime using only the ip^ boundary 
wave from all three edges, the lifetime is tte ~ UTE/iPx- 
One finds 



wtte 



%/3 



2 W 1 — (sin 0c/ sin Y 



4 V c* / 



sm 



cos^ 9q 



— sm 



— I cos 



2 a- 





•(66) 



The second line of this formula highlights the difference 
for TE polarization compared to TM. The factor in the 
brackets is some number less than 1; it contrasts the 
bracket which reduces to cos^ 9c in formula ( |6^ for the 
TM polarization lifetime when /i = /i'. The crucial dif- 
ference is the factor e/e' ^ 1 present here, compared 
to a similar factor ji' / ji ~ 1 for the TM lifetime for- 
mula. This is the more dominant factor, and it suggests 
that roughly speaking, the ratio of the lifetimes for the 
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FIG. 10: Estimated lifetimes for some low TM modes indi- 
cated by (m,n) pairs, versus the index ratio N. The lifetime 
is scaled by the speed of light c* in the cavity medium and 
the cavity edge size a. 



FIG. 11: Estimated lifetimes for low TE modes indicated 
by (m,n) pairs, versus the index ratio N. These lifetimes are 
larger than the corresponding TM lifetimes (Fig. |l^ when N 
is sufficiently larger than the cutoff value for that mode. 



two polarizations, which have the same (approximately 
DBC) boundary conditions and frequencies, is 



ttm 



- = (- 
e' Vn' 



(67) 



The result holds as long as the index mismatch is ad- 
equately large compared to the cutoff value needed to 
stabilize that mode by TIR. Otherwise, at smaller index 
mismatch, the TM lifetime can be longer than the TE 
lifetime. 

Some results for ttm are presented in Fig. |o|, showing 
lifetimes as functions of the index mismatch for some of 
the lowest modes. When scaled by the triangle size and 
light speed in the cavity, the lifetimes increase abruptly 
above the TIR confinement limits, eventually increasing 
at a slower rate. For the modes shown, dimensionlcss 
frequencies Lua/c* are typically numbers greater than 10 
(see Fig. with the values of tc* /a also of the order of 
10. Combining these rough results, the mode lifetimes in 
units of the mode periods T are similar to 



r 

T 



UJT 10 X 10 



27r 



27r 



16. 



(68) 



Assuming that the lifetime estimates have included the 
dominant loss mechanism in the cavity, this large result 
for r/T indicates that the original assumption of a res- 
onance mode weakly confined by TIR should be a valid 
concept. Obviously, this holds far enough above the TIR 
confinement limits only, keeping in mind the approxi- 
mate nature of the Dirichlct boundary conditions that 
were applied. 

Comparitive results for txe for the same mode indexes 
are shown in Fig. n^. Due to the presence of the extra 
factor of e/e' w N, these lifetimes increase much more 



rapidly than the TM lifetimes. As an example, for mode 
(3,5), the TE lifetime is about 6 times longer than the 
TM lifetime at N = 8. On the other hand, for the mode 
(1,3), which has a much larger TIR confinement limit, 
the TE lifetime is only about 1/3 longer than the TM 
lifetime at N = 8. Clearly, modes with (m, n) indexes 
nearly the same, as in the form (m, m -f 2) with large 
m, require smaller index mismatch for TIR confinement, 
and will more closely follow the lifetime ratio determined 
strongly by the index mismatch, Eq. (|6^). 

In recent experiments |Q triangular semiconductor cav- 
ities with edges ranging from 75 to 350 /im were used. 
Assuming an effective index of refraction around n « 4, 
with vacuum on the exterior, and using a « 100/Ltm, Eqs. 
( |63| ) and (|6^) give rough lower estimates ttm ~ 1-3 ps, 
and Tte ^ 20 ps. Of course, for practical purposes of 
maintaining a resonating mode, a large value of t/T is 
much more relevant, as discussed above. Based on these 
calculations, smaller cavities will have reduced lifetimes, 
but also shorter oscillation periods in the same ratio. For 
large index ratio, where sin 0c ^ 1, cos^c ~ 1, and us- 
ing the fact that to < n, lifetime expression (^2|) and 
frequency expression (^) produce the estimate. 



T V3 



(69) 



This ratio is independent of the cavity size or dielectric 
properties, increasing only with the squared mode quan- 
tum number n. The TE mode lifetime under these as- 
sumptions should be even larger, by the squared refrac- 
tive index ratio (n/n')^. 
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V. CONCLUSIONS 

The resonant modes of an optical cavity with an equi- 
lateral 2D cross-section have been examined in an ap- 
proximate manner, starting from the exact analytic so- 
lutions for Dirichlet boundary conditions. The six plane 
waves that make up each mode have been analyzed for 
the conditions necessary such that all are confined in the 
cavity by TIR. For 2D electromagnetics (fcz = 0), modes 
with larger ratios of quantum indexes m/n are most eas- 
ily confined; conversely, modes with m = never undergo 
TIR confinement. The presence of a nonzero longitudinal 
wavevector kz can be expected to improve the possibili- 
ties for TIR confinement. 



The critical index ratios for TIR confinement are the 
same for TM and TE polarizations. The different Fres- 
nel factors, however, produce different rates of energy lost 
from the cavity by the evanescent boundary waves. For 
the 2D problem, this was shown to lead to considerably 
longer lifetimes of the TE polarization, enhanced approx- 
imately by a factor of the squared index ratio (n/n')^ 
conpared to the TM lifetime. The differences in these 
lifetimes would be expected to imply stronger coupling 
of EM fields from outside to inside the cavity in the TM 
polarization; it suggests that stimulation and generation 
of the modes by an external light source should be more 
efficient for TM polarization. 
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